Rayleigh— Taylor turbulence in two dimensions 
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The first consistent phenomenological theory for two and three dimensional Rayleigh-Taylor (RT) 
turbulence has recently been presented by Chertkov [Phys. Rev. Lett. 91 115001 (2003)]. By 
means of direct numerical simulations we confirm the spatio/temporal prediction of the theory in 
two dimensions and explore the breakdown of the phenomenological description due to intermittency 
effects. We show that small-scale statistics of velocity and temperature follow Bolgiano-Obukhov 
scaling. At the level of global observables we show that the time-dependent Nusselt and Reynolds 
numbers scale as the square root of the Rayleigh number. These results point to the conclusion 
that Rayleigh-Taylor turbulence in two and three dimensions, thanks to the absence of boundaries, 
provides a natural physical realization of the Kraichnan scaling regime hitherto associated with the 
elusive "ultimate state of thermal convection" . 

PACS numbers: 47.27.-i 



The Rayleigh-Taylor (RT) instability is a well-known 
fluid-mixing mechanism occuring when a light fluid is 
accelerated into a heavy fluid. For a fluid in a gravi- 
tational field, such a mechanism was first discovered by 
Lord Rayleigh in the 1880s 1] and later applied to all 
accelerated fluids by Sir Geoffrey Taylor in 1950 
RT instability plays a crucial role in many field of sci- 
ence and technology. As an example, large-scale mixing 
in the ejecta of a supernova explosion can be explained 
as a combination of the Rayleigh-Taylor and Kelvin- 
Helmholtz instabilities RT instability also plays a 
crucial role in inertial confinement fusion as it finally 
causes fuel-pusher mixing that potentially quenches ther- 
monuclear ignition. Suppression of the RT instability is 
thus very crucial for the ultimate goal of inertial fusion 
energy. The final stage of RT instability necessarily leads 
to the so-called Rayleigh-Taylor turbulence, the main 
subject of the present Letter. Despite the long history of 
RT turbulence, a consistent phenomenological theory has 
been presented only very recently by Chertkov Q. Dif- 
ferent behaviors are expected for the 3D and the 2D case. 
About the former regime, the "5/3"-Kolmogorov scenario 
is predicted, while the Bolgiano picture 6] is expected 
for the 2D case. This Letter presents the first attempt to 
compare numerical results with such phenomenological 
theory. We show that: (i) low-order statistics of temper- 
ature and velocity follow Bolgiano scaling; (ii) there are 
strong corrections (intermittency) for higher-order tem- 
perature statistics; (Hi) the behavior of time-dependent 
global quantities such as the Nusselt and Reynolds num- 
ber as a function of Rayleigh number follows Kraichnan 
scaling. 

The equations ruling the fluid evolution in the 2D 
Boussinesq approximation are: 

dtT + v-VT = kAT , (1) 
dtuj + v-Vuj = vAuj - /3VT X g , (2) 



T being the temperature field, a; = V x t; the vorticity, 
g the gravitational acceleration, /3 the thermal expansion 
coefficient, k molecular diffusivity and v viscosity. 
At time t = 0, the system is at rest with the colder fluid 
placed above the hotter one. This amounts to assum- 
ing a step function for the initial tcmpertaure profile: 
T(0, a;) = — s(7n(z)8/2, 8 being the initial temperature 
jump. At sufficiently long times a mixing layer of width 
L[t) sets in, giving rise to a fully developed, nonstation- 
ary, turbulent zone, growing in time as L(f) ^ (i-e., 
with velocity uj^it) ^ t). If, on one hand, there is a 
general consensus on this quadratic law, which also has 
a simple physical meaning in terms of gravitational fall 
and rise of thermal plumes, on the other hand the value of 
the prefactor and its possible universality is still a much 
debated issue (see, e.g. Ref. .7J and references therein). 
The statistics of velocity and temperature fiuctuations 
inside the mixing zone is the realm of application of the 
pheomenological theory of Ref. 4]. Let us briefiy re- 
call the main predictions of this theory and some of its 
merits and intrinsic limitations. The cornerstone of the 
theory is the quasi-equilibrium picture where small scales 
adjust adiabaticafly as temperature and velocity fiuctu- 
ations decay in time. Upon assuming that the temper- 
ature behaves as a passive scalar, the analysis of two- 
dimensional Navier-Stokes turbulence leads to two sce- 
narios. While temperature variance flows to small scales 
at a constant flux, the velocity field either undergoes an 
inverse cascade with an inertial range characterized by a 
backward scale-independent energy flux or it develops a 
direct enstrophy cascade (for background information on 
two-dimensional turbulence see Ref. j3| for a theoretical 
introduction and Refs. for experimentally oriented 

reviews). Both possibilities actually turn out to be incon- 
sistent 01 ■ 

This apparent deadlock can be broken by rejecting 
the initial assumption that temperature behaves as a 
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passively transported quantity at all scales. Indeed, 
Chertkov suggests that buoyancy and nonlinear terms 
in Eq. (O must be in equilibrium. This is the essence 
of the Bolgiano regime [6j, and under the assumption 
that temperature fluctuations cascade to small scales at 
a constant rate one arrives to the Bolgiano scaling 
relations: 

s^T ^e(;^)'^' - (/35)^'/'^e4/v/5t-2/5 (3) 

S^v ^ ULit) [j^y^' - iPgef/'^ryH-'/' . (4) 

The above results constitute a set of mean field (i.e. di- 
mensional) predictions which need to be verified against 
numerical simulations and/or experiments. Our aim here 
is to shed some light on both the theory proposed by 
Chertkov and to expose the presence of intermittent phe- 
nomena (that could not be addressed within the phe- 
nomenological framework of Ref . Q ) by means of direct 
numerical simulations of equations , Q . 
The integration of both equations is performed by a 
standard 2/3-dealiased pseudospectral method on a dou- 
bly periodic domain of horizontal/vertical aspect ratio 
Lx/Lz = 1/4. The resolution is 128 x 4096 collocation 
points. Different aspect ratios (up to 1 : 8) and resolu- 
tions (up to 128 X 8192) did not show substantial modifi- 
cations on the results. In order to avoid possible inertial 
range contaminations, no hyperviscosity/hyperdiffusivity 
have been used. The time evolution is implemented by a 
standard second-order Runge-Kutta scheme. The inte- 
gration starts from an initial condition corresponding to a 
zero velocity field and to a step function for the tempera- 
ture. Given that the system is intrinsically nonstationary, 
averages to compute statistical observables are performed 
over different realizations (about 40 in the present study). 
The latter are produced by generating initial interfaces 
with sinusoidal waves of equal amplitude and random 
phases • Each realization is advanced in time until the 
mixing layer invades the 75% of the domain. The hor- 
izontally ensemble-averaged temperature field at three 
different instants is shown in Fig. ^ It is worth noticing 
the almost linear behavior of the averaged temperature 
within the mixed layer. This is a first clue suggesting 
a possible relation between RT turbulence and the 2D 
Boussinesq driven convection studied in Refs. 0, H^ . 
Further evidences will be given momentarily. In that 
particular instance of two-dimensional convection, tur- 
bulent fluctuations are driven by an external, linearly 
behaving with the elevation, temperature profile and the 
emer genc e of the Bolgiano regime clearly appears from 
data hj. We will argue that 2D RT turbulence corre- 
sponds to the case driven by a linear temperature profile 
with a mean gradient that adiabatically decreases in time 
as Q/L^t-^. 

The mixing layer growth rate is shown in Fig. [3 in terms 
of the growth-rate parameter a. Consistently with previ- 
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FIG. 1: The liorizontally ensemble-averaged temperature field 
at t/r — 2.6 (solid line), t/r — 3.9 (long-dashed line) and at 
t/r = 4.4 (short-dashed line). The dimensional group r = 
{Lz /{Ag)y^^ is a characteristic time scale of the flow. Here, 
A is the Atwood number, g the gravitational acceleration and 
Lz is the vertical size of the box. The value of Ag is 0.15. Note 
the almost linear behaviour of {T}z in the mixed layer. 



ous findings the mixing layer grows quadratically in time, 
i.e. a becomes almost constant in time reaching the 
value a ~ 0.12. The latter value is in agrement with the 
one found in Ref. (see the case x = 0.01). 

In order to quantitatively assess the presence of 
Bolgiano's regime for RT turbulence, let us focus on 
scaling behavior of stucture functions. For both velocity 
and temperature differences they are shown in Fig. |31 
On the one hand, second-order structure functions 
follow the dimensional predictions Q and |0J. On 
the other hand, moments of order 4 and 6 display 
intermittency corrections for the temperature, e. g. 
deviations from (jj)), while this is not the case for the 
velocity which shows a close-to-Gaussian probability 
density function for inertial range increments (not 
shown). The slopes of Fig. |21we have associated to Sj 
and are relative to the scaling exponents found in 
Ref. hJ. This is a further quantitative evidence in 
favour of the equivalence between RT turbulence and 
Boussinesq turbulence in two dimensions 11]. With the 
present statistics, moments of order higher than 6 are 
not accessible. We can further corroborate our claim 
on the possible equivalence between RT and driven 
2D Boussinesq turbulence by looking at the temporal 
behavior of structure functions. Dimensional predictions 
immediately follow from Eqs. Q and Q). The latter 
are well verified for all displayed orders for the velocity 
field (see Fig. For the temperature field, anomalous 
corrections start to appear at the fourth order and are 
of the form S'J(r) ~ e"(r/L(i))"/5(r/L(i))-'^" . If we 
assume (see Fig. OJ that the present RT model possesses 
the same spatial scaling exponents as those of the model 
presented in Ref. [l^, i.e. Sj{r) - r^-^ and S^ir) - r"-^ 
(and thus (T4 = 0.2 and erg = 0.5), we immediately get 
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FIG. 2: The time evolution of the layer-growth parame- 
ter a = [{l/Ag)]dL/d{t'^), L being the mixing layer width. 
Heavy line represents the ensemble-averaged quantity while 
thin lines refer to each individual realization thus giving an in- 
dication of the level of fluctuations. The horizontal line refers 
to the value of a found in Q and it is shown for compari- 
son. The mixing layer width L{t) is defined as the distance 
between z-levels at which J = ((T)^ - Tc)/{Th - Tc) = 1% 
and 99%, respectively. Tc and Th are the temperature of the 
cold and hot fluid, before the mixing process takes place. Two 
snapshots of the temperature fleld are shown a.t t/r = 2.6 (on 
the left) and at t/r = 4.4 (on the right). Dark (white) areas 
identify cold (hot) regions. 

a prediction for the exponents relative to the temporal 
behavior. We just have to remember that L{t) ^ 
to obtain the scaling relations ^4 ~ i^"-^ ~ t^^ '^ and 
~ L^"'^ ^ t~^'^. The latter are compatible with our 
results presented in Fig.^(b). 

We end up by discussing the behavior of turbulent heat 
flux, mean temperature gradient and root-mean-square 
velocity as a function of time. These quantities are cus- 
tomarily represented in adimensional variables by the 
Nusselt number Nu = [{vzT)L/{Ke)] + 1, the Rayleigh 
number Ra — g(3QL^ / {vk) and the Reynolds number 
Re = UrmsL/v. The question about what functional re- 
lations exist among these quantities is a long-debated is- 
sue in the context of three-dimensional Rayleigh-Benard 
turbulence (see, for example, Refs. fli. EE ES 113, Ell 
and citations therein). In 2D RT turbulence exact ex- 
pressions linking these adimensional numbers to temper- 
ature and kinetic energy input/dissipation rates can be 
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FIG. 3: The moments of longitudinal velocity differences 
SXif) (a) and temperature differences S'^{r) (b). We con- 
sider the isotropic contribution to the statistics, by averaging 
over all directions of separations r. In (a) dashed lines are the 
Bolgiano's dimensional predictions, SX ^ r^"/^ [gj. In (b), for 
n = 2 the dashed line is the Bolgiano's dimensional prediction: 
S2 ~ r^^^ @. For n = 4, 6 scaling exponents are anomalous, 
and their values are compatible with those of the 2D turbulent 
convection model of Refs. [TH (dashed lines). Moments 
are averaged over different L's ranging from L/L^ = 0.4 to 
L/L, = 0.6 (Fig.BJ. 



derived closely following Ref. [15|. For the RT case we 
have to deal with the additional complication of time- 
dependence yet compensated by the simplification orig- 
inating from the absence of boundary effects. These re- 
lations read dt{v'^/2) = vk^L^'^{Nu — l)Ra — and 
df{0^/2) = kQ'^L-\Nu - 1) - ee, being the depar- 
ture from the mean temperature profile. Since in 2D 
RT kinetic energy is transferred upscale we have a neg- 
ligible e„ and we can estimate a rate-of-change of ki- 
netic energy (PyQ^t from Eq. Q and L — PgQt^. 
For the temperature we have (6*^) ~ 8^ independent 
of time. Temperature performs a direct cascade and 
thus dissipation can be estimated as 5l0^5lv/L ^ 
Plugging those estimates into the exact relations yields 
Nu - {[ig&fK-^t^ and Ra ~ {(igQY{yn)-'^t^ , therefore 
Nu ^ Ra^/'^Pr^/'^. As for the Reynolds number, Eq. Q 
gives Re ^ {(igQfv-^t^ ^ Ra^/'^Pr~'^/'^ . It is worth re- 
marking that a similar analysis for the 3D RT case leads 
mutatis mutandis to the same scaling laws. These, when 
considered as a function of i?a, coincide with the results 
derived by Kraichnan for the pure bulk contribution to 
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FIG. 4: Evolution with L of (a) moments of velocity differ- 
ences SX and (b) temperature differences S^. The results 
are obtained by averaging over a fixed range of separations 
belonging to the inertial scales. We found a better scaling be- 
havior by displaying the behavior of structure functions as a 
function of the mixing layer width, L, rather than with respect 
to time. Those are connected by the relation L ^ t^. The 
dash-dotted lines correspond to the dimensional scaling for 
velocity and to the intermittency-corrected scaling for tem- 
perature (see text for details). 




FIG. 5: The behavior of the Nusselt number (bullets) and 
Reynolds number (circles) versus the Rayleigh number. 



3D Rayleigh-Benard turbulence. However, such Kraich- 
nan scaling regime for RB convection, also dubbed "the 
ultimate state of thermal convection" , has so far eluded 
both experimental and numerical confirmation [Toj . Ad- 
ditionally, it has been shown that it is not realizable in 
the analytically tractable case of Pr going to infinity . 
The reason may be traced back to the fundamental role 
played by the boundaries in establishing the turbulent 



heat transport in Rayleigh-Benard convection. Indeed, 
when boundaries are artificially removed as in the numer- 
ical simulations of Refs. the Kraichnan scaling is 
clearly observed. In this context, Rayleigh- Taylor turbu- 
lence provides a natural framework where heat transport 
takes place exclusively by bulk mechanisms and thus pro- 
vides a physically realizable example of the Kraichnan 
scaling regime, inviting further experimental and numer- 
ical effort in this direction. 
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